Efficient intrinsic spin relaxation in graphene due to flexural distortions 
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We propose an intrinsic spin scattering mechanism in graphene originated by the interplay of 
atomic spin-orbit interaction and the local curvature induced by flexural distortions of the atomic 
lattice. Starting from a multiorbital tight-binding Hamiltonian with spin-orbit coupling considered 
non-perturbatively, we derive an effective Hamiltonian for the spin scattering of the Dirac electrons 
due to flexural distortions. We compute the spin lifetime due to both flexural phonons and ripples 
and we find values in the 1-10 ns range at room temperature. The proposed mechanism dominates 
the spin relaxation in high mobility graphene samples and should also apply to other planar aromatic 
compounds. 

PACS numbers: 



Introduction. The electron spin lifetime in carbon 
materials is expected to be very long both because of 
the very large natural abundance of the spinless nuclear 
isotope ^^C and the small size of spin orbit coupling. 
In the case of flat graphene, the spin projection perpen- 
dicular to the plane is conserved, even in the presence 
of the intrinsic spin-orbit coupling. Thus, graphene was 
proposed as an optimal material to store quantum infor- 
mation in the spin of confined electrons [1] . Experiments 
PHZ] show that the spin lifetimes are much shorter than 
expected from these considerations, which is convention- 
ally attributed to the combination of two extrinsic fac- 
tors: the breaking of reflection symmetry due to coupling 
of graphene to a substrate 8J and/or a gate field, and the 
breaking of translational invariance due to impurities 19] . 

Here we propose a radically different scenario: efficient 
spin relaxation is intrinsic to graphene due to the inter- 
play between its unique mechanical and electronic prop- 
erties. We show that flexural distortions, unavoidable in 
two dimensional crystals [TU] in the form of either static 
ripples or phonons, induce spin scattering between the 
Dirac electrons, to linear order in the flexural field. This 
coupling differs from the spin-conserving second-order in- 
teraction between Dirac electrons and flexural distortions 
[TT] that has been proposed as an intrinsic limit to mo- 
bility in suspended high quality samples [12l [13] . 

The fact that curvature enhances spin-orbit scatter- 
ing has been discussed in the case of carbon nanotubes 
[14], graphene p^HlT] and graphene ribbons [18]. In 
the present Letter we derive a microscopic Hamiltonian 
that describes explicitly the spin-flip scattering of elec- 
tronic states of graphene due to both dynamic and static 
flexural distortions. We describe graphene by means of 
a multi-orbital atomistic description that naturally ac- 
counts for the two crucial ingredients of the proposed in- 
trinsic spin-phonon coupling: the intra-atomic spin orbit 



coupling (SOC) and the modulation of the inter-atomic 
integrals due to atomic displacements. Importantly, both 
the SOC and the flexural distortions couple the Dirac 
electrons to higher energy a bands, in the spin- flip and 
spin conserving channels respectively. Their combined 
action results in an effective spin-flip interaction for the 
Dirac electrons. The spin-flip scattering rates computed 
from the theory are shown to provide a very efficient in- 
trinsic spin relaxation mechanism in graphene at room 
temperature, without the need to invoke scattering with 
impurities or extrinsic SOC. We find that the spin lifte- 
time of electrons depends crucially on the long wave- 
length mechanical properties of the sample, determined 
by its coupling to the environment. 

Microscopic model. Our starting point is the tight- 
binding Hamiltonian H = Hsc + T~Lsoc for electrons 
with spin s moving in a lattice of atoms r, with atomic 
orbitals a. We write the hopping part as 

nsC= E H,.o'{r-^H^o,s'^r',o',s- (1) 
r,f' .o.o' ,s 

considering explicitly the dependence of the (spin con- 
serving) inter-atomic matrix elements on the positions of 
the atoms. The intra-atomic SOC reads 

-Hsoc^ A(foa|L(f) • S'lro's')*!,,,,*.-^',.' (2) 

r,o,o',s,s' 

where S is the spin operator, L{f) is the orbital angular 
momentum operator acting upon the atomic orbitals of 
site f and A is the spin-orbit coupling parameter. 

Deviations from the ideal graphene lattice affect its 
electronic properties via modifications of the transfer in- 
tegrals in Eq. ([T]) . Their dependence on the inter-atomic 
distance, for example, gives rise to an electron-phonon 
interaction analogous to that of conducting polymers 
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[T^ and leads to the appearence of effective gauge fields 
[TBI [20] . In the present model, we consider instead the 
coupling with flexural distortions arising from the angu- 
lar dependence of the interatomic Hamiltonian. We de- 
scribe corrugations away from perfectly flat graphene in 
the form r ~ rQ + h{f)z, where ro is a vector of the honey- 
comb lattice and h{r) is the displacement of atom r per- 
pendicular to the graphene sheet. We expand the inter- 
atomic Hamiltonian matrix to lowest order in the flexural 
field h{r) and rewrite the Hamiltonian as H = Hq + V, 
where "Ho now describes ideally flat graphene including 
the weak intra-atomic SOC perturbation, and 

r,r' ,o,o' .s 

(3) 

is the spin-conserving coupling between electrons and 
corrugations. 

Electron- flexural phonon scattering. It is now conve- 
nient to recast Eq. ([S]) in terms of the eigenstates of the 
Hamiltonian T-Iq for flat graphene. In this context, atoms 
are specified by their unit cell vector, R, and their posi- 
tion Tq, inside the cell (sublattice index a = A, B). The 
creation operators for Bloch states are related to atomic 
orbitals through: 



1 



N J-^ 



ik-R 



(4) 



where k is the wave-vector, the coefficients j^{a^o, s) 
are obtained from the diagonalization of the Bloch ma- 
trix and v is an index that labels the resulting bands 
(with mixed spin and angular momentum). Similarly, 
we expand the flexural fleld on each sublattice in its 
Fourier components, ha{R) = -^J2q^'^^^^^a{q)- Af- 
ter a lengthy but straightforward calculation, we can ex- 
press Eq. ([3]) as a term causing scattering between crystal 
states with different momentum and band index: 



(5) 



^ ^ 1 d 

,,'(k,k') = V TrHo.o'ira-ra' - R) 'X 

y/N ^ OZ 
R,a^a' ,0,0' ^(T 

X F£,{k,k')Cl^^{a,o,a)C^,j:,{a\o',a) (6) 



The coupling is linear in the flexural phonon field, 
through the form factor F^^,{k, k') = ha{k — k')e'^^ — 
ha'{k — k')e^'^'^. This should be contrasted with 
the electron-flexural phonon coupling usually considered 
within the tt subspace [TlllT3l [20]. that is quadratic in the 
fleld because of the quadratic dependence of interactomic 
distances on h. Note that i?oo' is short ranged within our 
tight-binding description, which limits R to the 4 vectors 



connecting neighboring cells (inter-cell coupling), as il- 
lustrated in Fig. [iji, plus the null vector (intra-cell cou- 
pling). 
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FIG. 1: (a) Sketch of a local corrugation of the graphene 
layer. The shaded area is the unit cell. Green arrows indicate 
the inter-atomic contributions to the electron-flexural phonon 
coupling of Eq. Q. (b) The band structure resulting from 
our Slater-Koster parametrization (c) Kinematics of the scat- 
tering process around the Dirac point K in the low phonon 
frequency limit. 



Slater-Koster parametrization. Equations (5[6) pro- 
vide a general recipe to compute the coupling of electrons 
of a generic tight-binding Hamiltonian to a flexural field. 
We now show that because we have included the SOC in 
the reference Hamiltonian Tio, the perturbation Eq. ([5| 
is able to induce a direct coupling between states with op- 
posite spin. Following previous work [T31 [121 [IZl [T51 [55] 
we consider a subset of 4 valence orbitals of the Car- 
bon atom, namely o = s,Px,Py,Pz, and adopt a Slater- 
Koster (SK) f2T] parametrization for the tight-binding 
Hamiltonian Eq. ^ [53 . In this framework, the inter- 
atomic matrix elements Hq o' connecting the atom a at 
ra with atom a' at fa' + R can be expressed in terms 
of 4 parameters, Vss, Vsp, Vc-k, F^tt, describing inter- 
orbital overlaps in the s,p basis [iHl [21], and the three 
director cosines l,m,n of the interatomic bond vectors 
p — fa — fa' ~ R, defined hy p = p (Ix -\- rwy -\- nz). Set- 
ting Ux = VppTT + x'^Van and Va^r = (Vppa - Vppjr) WC Cau 

write the SK matrix in compact form as 



H{P) 



ssa 



\ -nVs 



spa 



^^spa 

Ui 

ImVcrn 

InVcTTT 



mVspa 
ImVaj, 



nVspa 

InV^T, 
mnVa-TT 



(7) 



The unperturbed crystal states for flat graphene are 
described by H{p) with n = (all bonds within the x, y 
plane). The resulting band structure is shown in Fig. 
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[T|d. From Eq. the electron-flexural distortion cou- 
pUng is determined by dzH{p) = {l/d)dnH{p), where d 
is the equihbrium C-C distance. Direct inspection of Eq. 
Q shows that H{p){n = 0) does not couple the Pz and 
s,Px,Py sectors, while dnH{p) does. As a result, the tt 
and a bands of flat graphene are not mixed in the absence 
of SOC, unless they scatter with flexural distortions. In 
the presence of SOC, however, the tt states with spin s 
mix with the a states with opposite spin already within 
the reference "Hq for A^^t graphene. Close to the Dirac 
points, where the tt and a bands are separated in en- 
ergy by a gap Ea-n , the spin vr — cr mixing is proportional 
to A/i?CTir- Since this correction is small, the low energy 
Dirac bands of "Ho can still be labeled according to their 
dominant spin character, that we denote as and JJ-. 

The above derivation shows that there are two pertur- 
bations that couple tt and cr states, the spin-conserving 
coupling to the flexural field, and the spin-flip SOC. Their 
combination is able to yield a spin-flip channel within 
the low energy tt bands, that is linear in both the flexu- 
ral deformation and in the atomic spin-orbit coupling A, 
in contrast with the perturbation that applies to flat 
graphene |15j . 

Spin-flip rates. We now apply the microscopic theory 
developed above to compute the spin relaxation rate for 
electrons close to the Dirac points. Anticipating that the 
dominant contributions to spin-flip scattering arise from 
the long- wavelength, low energy flexural modes, we con- 
sider the lower flexural branch for which ojq cx (see be- 
low) , discarding the higher energy modes that involve out 
of phase vibrations of the two sublattices. The flexural 
fleld is factored out from Eq. ([5| by setting /ia('z) = 
and hsiq) ^ e''^<^''~^-^^h^, which yields: 



(8) 



with Vy ,y'(fc + (f, fc) — -^M^_^t{k,q). The standard 
form for the phonon spin-flip interaction in second 
quantization is readily obtained by substituting — 



-(at- -I- flg-) into Eq. 



(8), with Mc the Carbon 



mass. 

The spin relaxation rate can now be calculated from 
Eq. ([s]) via the Fermi golden rule. Because the dispersion 
of the flexural modes is much weaker than the electronic 
dispersion, we can safely neglect the phonon frequency 
in the energy conservation. The relaxation rate for an 
electron with momentum k in the band is then obtained 
by summing over both phonon absorption and emission 
processes and over all possible final states in the JJ- band, 
which yields: 



r,-: = 



27r 



d^q 



^ - h J j2^2\Ml^,^i^M\hf)5{E^,^^- E^). (9) 

This, together with the explicit expressions for the spin- 
flip matrix elements in Eq. ([6]), constitutes the main 



result of this work. From Eq. ^ it is clear that once 
that the specific form of the electron-flexural phonon cou- 
pling, M||._j|(A;,q), is known, the behavior of the spin 
relaxation rate is fully determined by by the statisti- 
cal fluctuations of the flexural field, (/i?.). Interestingly, 
the above expression describes on an equal footing both 
low-frequency flexural modes (that arise in free-standing 
or weakly bound graphene or graphite) as well as static 
ripples (relevant to graphene deposited on a substrate). 
The proposed spin relaxation mechanism therefore ap- 
plies without distinction to both physical situations. 

For actual calculations we approximate the tt band en- 
ergies as Ej: — -iihvF\k\, which is valid except for a neg- 
ligible interval around the Dirac point, where the SOC 
opens a gap of the order of few ^eV . Energy conservation 
implies \/k'^ + q^ + 2kq cos 9 = Vk^. This fixes the rela- 
tive angle between k and q to cos = (see Figjlj;) , 
which allows to perform the angular integration in Eq. 
^, yielding: 
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2k 



dq 



hvpy/l - {q/2kf 



(10) 



In principle, the spatial derivative in the electron- 
phonon coupling elements of Eq. ([6| should translate 
into a linear scaling of the matrix elements with ex- 
changed momentum q. This result, however, only holds 
at the Dirac points, i.e. for A: = 0. From our numerics, we 
find instead that in the relevant case of small but finite k 
the matrix elements evaluated on the energy-conserving 
surface satisfy the non-trivial result: 

\M^,y^{k,q)\^ = \\kdf[-l^)\ (11) 

independent on the SK parameters. 

Fluctuations of the flexural field. We now evaluate 
the flexural fluctuations (/i^) in different cases of physical 
interest. We start with the expression for free-standing 
graphene at thermal equilibrium: 



[l + 2nB{ujg)] 



(12) 



where nB{uOq) is the thermal population of mode q and 
the second equality holds when ksT ^ hujq. For purely 
harmonic flexural modes, for which ujg ~ Dq^ , the fluctu- 
ation (hi) diverges as q~* for small q. This divergence is 
renormalized due to the interaction with other phonons 
(i.e. by anharmonic effects) and can be further cut 
off by the presence of strain [13; or pinning to a substrate 
|25| . The resulting dispersion can be written in general 
as ujq = D\/q^ + q^~^^q2 so that, in the long wavelength 
limit, 



(hi) 



q-^-'iqi 



(13) 
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where ry and qc depend on the physical mechanism of 
renormalization. Specifically, substrate pinning opens a 
gap in the phonon spectrum [25|. corresponding to 77 = 4; 
strain makes the dispersion linear at long wavelengths 
[I3j (r; = 2); anharmonic effects yield 77 = 0.82 P] , 
Based on the asymptotic long- wavelength behavior Eq. 



n(10'^ cm"^) 
0.1 12 4 



( 13 1 we obtain the following general expression for the 



spin relaxation rate in terms of the scaling parameter 77, 
momentum k, and cutoff valid for k < q^. 



Tk = 



A2 {hl/^ 



El. d? Sirh 



(14) 



where = (V7r/2)r[(l + 77)/2]/r[l + 77/2] is a constant 
of order 1 and {hl/^)/cP = kBT{cP /McD) is a dimen- 
sionless measure of the C-C bond angle fluctuations [from 
Eq. (12), evaluated at g = 1/d S> gc]- The spin relax- 



ation rate increases linearly with temperature, following 
the thermal population of flexural phonons. A weaker 
temperature dependence arises when anharmonic effects 
dominate, because the anharmonic cutoff is itself tem- 
perature dependent, qc oc VT [lOl [24]. 

Interestingly, substituting the Fermi wavevector kp = 
y/wn into the above expression results in a density de- 



pencence cx n 



(which is valid for Ep > ksT). 



Therefore, the spin relaxation rate can either be an in- 
creasing or decreasing function of doping depending on 
the value of 77, which is ultimately determined by the 
mechanical environment of graphene. At high density 
{kp ^ Qc), instead, the cut-off becomes negligible and 
the integral in Eq. ( 10 1 converges to the result for ideal 
graphene [Eq. (14) with ?7 = 0] regardless of the value 
of rj. In view of these results, and because the proposed 
flexural spin relaxation mechanism is essentially indepen- 
dent from the spin- conserving scattering that determines 
the electrical mobility, some care should be taken when 
identifying the spin relaxation mechanism (Elliot- Yafet, 
or Dyakonov-Perel) from the relative scaling of the spin 
and momentum relaxation rates with carrier density, as 
is customarily done [51 [71 [THl [31] • 

Results and discussion. The spin lifetime = 1 /Tkp 



obtained from Eq. ( 10 ) for electrons at the Fermi level 
is plotted as a function of electron density in Fig. [2] 
We take D = 4.6 • IQ-'^m^s-^ [T^, A = 8 meV [25], 
T = 300K and vp = 1.16 • lO^m/s from our SK band 
structure. Different curves correspond to different val- 
ues of the scaling exponent rj, i.e. to physically different 
mechanical environments for graphene. Two representa- 
tive values of the cutoff momentum, q^. = 0.01 and 
qc = 0.1 A~'^ are considered, covering the large spread of 
qc values available in the literature. In both panels, the 
result for ideal graphene in the harmonic approximation 
is shown for reference (black) as it provides an absolute 
lower bound to the actual lifetime. 

In the density range relevant to current graphene-based 
devices, the predicted spin lifetime is typically in the 
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FIG. 2: Room temperature spin lifetime calculated from the 
theory for electrons at the Fermi energy, for 2 different values 
of cutoff momentum: (a) = 0.01 and (b) qc = 0.1 

(right) and different long-wavelength scaling laws: ideal 
graphene (black, rj = 0); including anharmonic effects (red, 
rj — 0.82); including strain (orange, rj — 2) and substrate 
pinning (gray, 77 = 4). The black dashed line is for 77 = 1, 
which is representative for substrate roughness. 



range of 1-10 ns at room temperature. Because the 
spin relaxation is dominated by the low energy fluctu- 
ations of the membrane, the shortest spin lifetimes are 
obtained in suspended unstrained graphene (red curve), 
i.e. when external mechanical influences are minimized 
and the mobility is possibly largest. This should also ap- 
ply to the case of graphene on boron nitride [27], where 
the substrate is atomically flat and effectively decou- 
pled from the graphene layer. All in all, considering 
the electronic mobilities currently reachable in suspended 
graphene /i = 2 • lO^crn^ /Vs (6 • IQ'^cin? /V s in the case 
of BN substrates [27]) our theoretical predictions for the 
spin lifetime set an intrinsic limit to the spin relaxation 
length in excess of Lg « lO/zm at room temperature. 

Longer spin lifetimes could be achieved by suppress- 
ing the fluctuations of the graphene membrane, by an 
applied strain (orange) [13] or substrate pinning (gray) 
(similarly, the proposed mechanism would be clearly in- 
effective in epitaxially grown graphene due to the strong 
interlayer binding) . We note however that while quench- 
ing dynamical fluctuations, the coupling to a substrate 
may introduce static corrugations, which are also detri- 
mental to the spin lifetime. This could be at the origin 
of the spin lifetimes measured in graphene on Si02 in 
Ref . [6j . Recent measurements of ripples of graphene on 
Si02 are indeed compatible with 77 < 1 for g > q* = 0.02 
A"-'^ [28], which is ascribed to the presence of substrate 
roughness. At densities such that kp > q* , the spin re- 



laxation from such ripples is again described by Eq. ( 14 ), 



with now a prefactor {hf^j)/d'^ determined by the sub- 



strate morphology rather than by Eq. (12), that could 
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in principle be extracted from the experiment. 

Concluding remarks. In summary, we have shown 
that corrugations, that are ubiquitous in graphene in the 
form of dynamical flexural phonons or static ripples, en- 
able a direct spin-flip mechanism to linear order in both 
the flexural field and in the spin-orbit coupling. This 
mechanism provides an efficient source of spin relaxation, 
setting an intrinsic limit for the spin lifetime in high mo- 
bility samples. Such limit is however non-universal, as 
its precise value depends on graphene's mechanical envi- 
ronment, that determines the long- wavelength behavior 
of the flexural field. At room temperature, intrinsic spin 
lifetimes in the nanosecond range are expected in a very 
wide range of situations. We note that an analogous 
mechanism could be at the origin of the spin relaxation 
recently reported in small organic molecules |32j . with 
an enhancement of the spin-orbit coupling caused by the 
low-frequency molecular bending modes. 

Finally, whereas the existence of a rather small upper 
limit for spin-lifetimes in graphene might present an ob- 
stacle for certain applications like spin transistors, the 
intrinsic spin-lattice coupling could open the way for hy- 
brid devices, where a confined vibrational phonon could 
be coupled resonantly coupled to the spin-flip transitions 
of Zeeman split of confined Dirac electrons. Microwave 
pumping of such system could result in a maser behavior 
of the phonon mode [55] . 
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